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Abstract 



AMS class.: 15-04; 15A06; 15A39. 
Keywords and phrases; tropical linear system, algorithm. 
An algorithm to give an explicit description of all the solutions to any trop- 
ical linear system AQx = BQxis presented. The given system is converted 
into a finite (rather small) number p of pairs {S, T) of classical linear systems: 
#y^ ] a system S of equations and a system T of inequalities. The notion, intro- 

^ ■ duced here, that makes p small, is called compatibility. The particular feature 

CO I of both S and T is that each item (equation or inequality) is bivariate, i.e., it 

O^ ■ involves exactly two variables; one variable with coefficient 1, and the other 

one with —1. 5 is solved by Gaussian elimination. We explain how to solve 
T by a method similar to Gaussian elimination. To achieve this, we introduce 
the notion of sub-special matrix. The procedure applied to T is, therefore, 
called sub-specialization. 



1 Introduction 



H I Consider the set M U {— cx)}, denoted T for short, endowed with tropical addition 

® and tropical multiplication 0, where these operations are defined as follows: 

a ® 6 = max{a, b}, aQ b = a + b, 

for a, 6 G T. Here, — oo is the neutral element for tropical addition and is the 
neutral element for tropical multiplication. Notice that a © a = a, for all a, i.e., 
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tropical addition is idempotent. Notice also that a has no inverse with respect to 
©. We will write © or max, (resp. or +) at our convenience. In this paper we 
will use the adjective classical as opposed to tropical. Most definitions in tropical 
mathematics just mimic the classical ones. Very often, working with (T, ©, 0) leads 
to working with min, which will be denoted ©'. 

Given matrices A,B^ A^mxn(T), we want to describe all x £ T"- such that 
AQ X = B Q X. This means 



max{aij + Xj : 1 < j < n} = max{6jj + Xj : 1 < j < n}, i = 1,2, 



,m. 



Of course, Xj = — oo, all j = 1, 2, . . . , n is a solution (the trivial solution). Using 
the notions of winning pairs, compatibility and win sequence (introduced in this 
paper; see. definitions [2j |7] and [TT]) . the given problem is reduced to solving a 
finite, rather small, number of pairs (5, T) of classical linear systems: a system S 
of equations and a system T of inequalities, each item (equation or inequality) being 
bivariate, i.e., it involves exactly two variables, one with coefficient 1, and another 
with coefficient —1. Of course, S can be easily solved by Gaussian elimination. 
On the other hand, the fact that T consists of bivariate inequalities allows us to 
suggest a Gaussian-like procedure to solve system T. More precisely, by certain 
row operations (see p. [T4l) . we do not triangulate a coefficient matrix for T, but 
nearly so. What we do is to transform such a matrix into two matrices, one of which 
is a sub-special matrix (see definition [T5l) . From these two matrices the solution 
set to the T is directly read. We call this procedure sub-specialization, (see remark 
[T8] ). Notice that we need not use the simplex algorithm or other well-known ones 
to solve T. 

Compatibility of winning pairs (see definition |7]) turns out a very handy nec- 
essary condition to work with. Indeed, two given winning pairs yield some bi- 
variate linear inequalities. Imagine that two of these are xi — X2 + a < and 
—xi + X2 — b <^, for some a, 6 € T. Now the compatibility of the winning pairs 
guarantees that these inequalities hold simultaneously, i.e., they can be concatenated 
into xi + a < X2 < xi + &, so that a < b. Therefore, by requiring compatibility 
of winning pairs, all we are doing is ruling out, at the very beginning, wasting out 
time with systems of inequalities which, for sure, the only solution is trivial. The 

problem AQ x = B Q x has been addressed before. Indeed, in lH, it is proved that 
the solution set can be finitely generated. In Q, a strongly polynomial algorithm is 
found which either finds a solution or tells us that no solution exists. In |lj sec. 3.5, 
all the solutions are computed by a technique called symmetrization and resolution 
of balances. In lH, generators for the solution set are computed. The idea of finding 
(a minimal family of) generators is pursued in |[T2l [141 [151 for the equivalent prob- 
lem AQ x < B Q x. An iterative method in presented in jQ for another equivalent 
problem, namely A Q x = B Q y, where x and y are unknown. Also, there is a 
technique in [TQ to solve the problem AQx(Ba = BQx®b, relying on a recursive 



formulation of the closure operator (also called Kleene star operator) on matrices. 
In [91 ch. 4, the closely related problem A Q x ® b = x (similar to the classical 
Jacobi iterative method) is solved using Kleene stars. 

We present the solution setto AQ x = B Q x as a. finite union of sets, which ai^e 
(attending to their presentation) obviously convex. This, of course, agrees with the 
cellular decomposition in fE\. Convexity issues are also studied in |[T3l . The recent 
paper ||T4]| addresses the problem Aq x < B Q x. The drawback of the algorithm 
presented in lfT4l is, in our opinion, that this algorithm calls for the calculation of 
certain Kleene stars, and this is not an easy task. 

We do not describe the solution set by generators. We do not use Kleene stars. 
Instead, our method provides a description of each convex piece of the solution set 
by parameters and bivariate linear- inequality relations among them, in a reduced 
form. 

Let m,n,s G N be given. The following are kindred problems in tropical linear 
algebra: 

• PI: AQx = 0, 

• P2: AQx = b, 

• P3: AQx <b, 

• PA: AQx = B Qx, 

• P5: AQx < BQx, 

• P6:CQx = DQy, 

• P7: AQx®a = BQx®b. 

Here the data are matrices A,Be MmxnO^), C € MsxnO^), D G Msxmi^) and 
vectors a, b over T, and the jth problem is computing all vectors x S T", y E T'", 
such that Pj holds. By deciding the jth problem we mean either finding one 
solution or declaring that the problem has no (non-trivial) solution. Some references 
since 1984 are |[T]|4l|6l|7J[12j[T5j[l6l. Eai^lier books and papers can be found there. 

Of course, x = — oo, y = — oo are solutions to P3, P4, P5 and P6. These are 
the trivial solutions. 

Only if the vector b is real, problem P2 reduces to PI. More generally, one 
must realize that, contrary to classical linear algebra, problems P7 and P4 do not 
reduce to problem P2 or PI, because there are no inverses for tropical addition, 
so there is no tropical analogue for the matrix —B. Nevertheless, there are well- 
known connections among these problems, i.e., being able to solve some of them is 
equivalent to being able to solve some other. 

We need some notations: 



• For c, d G T, c ©' d means min{c, d} and cQ' d means c + d. 

• For c,d ^ T", c ©' d'^ means min{ci + di, C2 + ^2, • • • , c„ + dn}- 

• If ^4 = {ttij) G Mmxn{^) then A* = {—aji) is the conjugate matrix. 
The relationship among these problems is as follows: 

• Deciding P3 is possible, if A is real. 

Indeed, x* = ^* 0' 6 is a solution (called principal solution) and x < x* 
if and only if x is a solution; see ||6l, p. 31; in HI this process is called 
residuation. 

• Deciding P3 helps with deciding P2, if A is real. 

Indeed, P2 might be incompatible but, if it has a solution, then x* is the 
greatest one; see 10, p. 31. 

• Deciding P6 implies deciding P2. 

Given A and b, we decide Aq x = I Q y. For each pair of solutions x, y, if 
any, we set y = b,if possible. 

• Deciding PA is equivalent to deciding P6. 

Suppose X is a solution to PA and write AQ x = y. Concatenating matrices, 
write C = ^ e 7W2mxn(T), D = { G X2mxn(T), whcrc / is the 
tropical identity matrix, so that C Q x = D Q y. Therefore, if we can decide 
P6, then we can decide PA. 



Suppose now x, y are solutions to P6 and write z = ^ , A = [C,—oo], 

B = [—oo, D] so that AQ z = B Q z. Therefore, if we can decide PA, then 
we can decide P6. 

• PA and P5 are equivalent. 

AQx = BQx is equivalent tO740x<i?0x and ^4 x > i? x. On the 
other hand, 740x<i?0xis equivalent to (^4 ® B) x = i? x. 

• Deciding PA implies deciding P7. 

We introduce a new scalar variable z and write AQxQaQz = BQxQbQz. 
Concatenating matrices, write t = [^], C = [A,a\, D = [B,b] so that 
C Qt = DQt. After solving PA, set i„+i = z = 0. 



2 The problem 

Given matrices A,B^ A4mxn(T), we want to describe all non-trivial x E T" such 
that 

AQx = BQx. (1) 

Notations: 

• A = {aij), B = (bij), with 





^^J = 


1 aij if aij >bij, u - \ ^»i ^^ "»J - ^«J' 
1— oo otherwise, 1 — oo otherwise. 


• M 


= A®B = 


= A © B = {rriij ) is the maximum matrix. 


Let 







A0x = B03; = M0x. (2) 

Notice that ^ is equivalent to (|2]) and © is simpler than ([T]) because it involves 
fewer real coefficients. Thus, we will assume that A = X and i? = B, {assumption 
i) in the following. 
More notations: 

• [n] = {1, . . . , n}, for n S N. 

• For any c G T, x = c S T" means Xj = c, for all j S [n]. 

• For convenience, we extend T to T, by adding +oo. Eventually we will 
remove +oo. 



If /c G [m] and j, / G [n], then 
dif(M;j,Ofc = 



rnkj-ruki, ifniki^-oo, 

+00 if mytj 7^ — oo and rriki = — oo, 

undetermined, otherwise. 



The undetermined case will never appear in the following. 

• = {j G [n] : Xj = — cx)}. 

We will deal with bivariate equalities and inequalities, linear on the Xj's with 
coefficients in T. Tautological linear equalities or inequalities will be always re- 
moved {assumption 2). A non-tautological linear equality not involving ±00 will 
be reduced to an equivalent equation of the form a = 0, for some a, by the usual 
algebraic rules. Of course, —a = is also possible. A non-tautological linear in- 
equality not involving ztcxo will be reduced to an equivalent inequality of the form 



a < 0, for some a. These will be called normal forms. Notice that normal forms 
have real coefficients. 

What will we do with certain normal forms (/>, if we know that j € Q, i.e., 
Xj = — oo? We will remove and enlarge as follows: 

• If (j)is Xj — Xk + a < 0, for some k ^ j and a G M, then remove (/>. 

• li 4) is Xk — Xj + a < 0, for some k ^ j and a G M, then remove (j) and set 
Xk = — oo, i.e., enlai^ge Q, with k. 



Remark 1. 



1. X 



-oo satisfies (|7]). This is the trivial solution. 



2. Ifiow{A,i) > row{B,i) or row(A,i) < iow{B,i) for some i G [m], then 
X = —oo is the only solution to (|7]). 

3. If row{A,i) = row{B,i) for some i G [m], then these two rows can be 
removed, so that m can be decreased to m — 1. 

4. Ifcol{A,j) = col{B,j) = —oo for some j G [n], then no restriction is 
imposed on Xj. Then these two columns and Xj can be removed, so that n 
decreases to n — 1. 

We will assume that row(^, i) ^ row(S, i) {assumption 3), row(^, i) ^ vow{B, i) 
(assumption 4) and row(^, i) ^z^ row{B,i) (assumption 5), for all i G [m], and 
col{A,j) = col{B,j) = — oo (assumption 6), for no j G [n], in the following. 

The sets in the next definition are denoted /. J.K.L in HI. 
Definition 2. For each i G [m], let 

1. WA{i) = {j : aij > b,j}, WB{i) = {j : Uij < hj}. 



2. E{i) = {j : aij = kj / -oo}, F{i) = {j : m 



V 



-oo}. 



3. win(i) = {WA(i) x WB{i)) U {E{i) x E{i)) C [n] x [n]. Each element of 
win(i) is called a winning pair. 

For each i G [m], WA{i) U WB{i) U E{i) U E{i) = [n] is a disjoint union. By 
assumption 3, E{i)\JE{i) / [n]. By assumptions 4 and 5, WA{i) / [n] / WB{i). 
Moreover, ^^^iF{h) = 0, by assumption 6. 

Example 3. Given 
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-1 8 


6 7 


5 1 
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1 2 



we get 

M 



Then WA{1) = {1,2,3}, WB{1) = {A}, WA{2) = {1,2}, WB{2) = {3,4}, 
WB{2,) = {A}andE{:i) = {1,2,3}. Therefore, win(l) = {(1,4), (2,4), (3, 4)}, 
wm(2) = {(1, 3), (1, 4), (2, 3), (2, 4)}, wm(3) = {(1, 1), (2, 2), (3, 3)}. 

Remark 4. Suppose that win(i) = 0. This is equivalent to either row(A, i) = — oo 
or j:ow{B,i) = — oo, but not both, by assumption 3. Say row(^,i) = — oo and 
bij 7^ —oo for some j € [n]. If A Qx = BQx = y, then yi = — oo and Xj = — oo. 
With this information in mind, we can remove the i—th rows, so that m decreases to 
771 — 1. Therefore, we will assume that row(j4, i) ^ — oo and iow{B, i) ^ —oo, for 
all i € [m] (assumption 7). Then win(i) ^ 0, for all i G [m]. 

Non-trivial solutions to ([T|l arise from winning pairs. Let us see how. Recall 
that M = A® B = {rriij) G Mmy,n{^)- M might not be real (see example]!]). 

Definition 5. Consider i G [m] and I G win(i). Let x G T", y G T™ be any 
vectors satisfying AQx = BQx = y(in particular, 

row(A, i) Q X = row{B,i) Q x = row(M, i) Q x = yi). 

We say that the solution x to O arises from / if 

ruij + Xj <yi, (3) 

for all j G [n] \ F{i), with equality for all j G |/|. 

Notice that j G F{i) if and only if ruij = — oo and, in such a case, the inequality 
© is tautological. Otherwise, rriij is real. 

Suppose that x arises from I and write |/| = {ii,i2}- Then, equality for all 
j G |/| means 

ITT'ii^ ~r Xi_^ — yi — ^^ii2 "•" '^12 ' 

whence we obtain one bivariate linear equation 

Xi^ = dif (M; ii , i2)i + Xi^ . (4) 

Notice that dif (M; 1^,^2)1 is real. In addition, ([3]) amounts to, at most, 2n — 2 
bivariate linear inequalities in the Xj 's. Notice that (|4| is tautological, when i^ = 12- 

Example 3. (Continued) Take i = 1, / = (1,4) G win(l) and suppose that a 
solution X to ^ arises from I. This means 

3 + xi = 8 + X4 



(one bivariate equation) and 

7 + X2<yi, -l + X3<yi, 

where 3 + xi = yi- Replacing yi by its value, we obtain two additional bivariate 
inequalities. Altogether, ((21) becomes, in normal form, 

xi — X4 — 5 = 0, (5) 

-xi + X2 + 4 < 0, (6) 

-xi + X3 - 4 < 0. (7) 

Remark 6. Suppose i,k (^ [m], i < k, I ^ win(i), K € win(A;). Assume that the 
solution X arises from I and from K. Then for all i € |/| and k £ \K\, we have 

mii + Xj^ = yi, (8) 

mik + Xk<yi, (9) 

mkk + Xk = yk, (10) 

mki + Xi_<yk. (11) 

Adding up, 

mik + mki + X£ + Xfc < rriii + mkk + Xj^ + Xk = yi + yk, 

whence 

rriik + ruki < mu + rukk- (12) 

In other words, the value of the 2x2 tropical minor of M, denoted M{i, k; i, k). 



rriii niik 
■mki rukk 



= max{mii + lUkk, rriik + mki}, (13) 

trop 



is attained at the main diagonal. One more way to write this over T is 

d[f{M;i,k)k < d[i{M;i,k)i. (14) 

This remark leads to the following key definition. 

Definition 7. Consider i,k G [m], i < k, I £ wm{i),K G wm{k). We say that 
K is compatible with I if the value of f UH ) is attained at the main diagonal, for all 

i £ \I\ and all k & \K\. Equivalently, if di4l) holds in T, for all i € |/j and all 

k£ \K\. 



Example 3. (Continued) In our running example, take i = 1, A; = 2, / = (1, 4) and 
K = (1,3). Then K is compatible with I, since each tropical minor 
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14, 



trop 



-1 



1 



13 



trop 



attains its value at the main diagonal. The inequalities di4D are < 0, 1 < 4, 
— 5 < 5 and — 4 < 9 in this case. However, K = (1,4) is not compatible with I, 
since the tropical minor 

3 5 



6 1 



14 



trop 



does not attain its value at the main diagonal. 

Remark 8. Suppose i,k (^ [m], i < k, I € wm{i),K € win(A;), K compatible 
with I. If i & \I\ and k G \K\ are fixed, then dii{M,i,k) is decreasing on the 
subscripts, by compatibility. Therefore 



[dif(M;i,Ai)fc,dif(M;i,Ai),] 



(15) 



is a non-empty closed interval, denoted mt{M;i,k;i,k). It may degenerate to 

[— oo, a], [a, +oo], [— oo, +oo] or [a, a] for some real a. 



amp] 


e 9. Given 




















A = 


1 

. a2i 


— oo — oo 
a22 


5 


B = 


— oo 

&21 


1 

^22 


— OO 




get 






















M = 


rr 


1 

^21 


1 - 

m22 


OO 




5 







some 021,022,^21, ^22, ?7i2i and 11122 S T. Take i = 1, I = (1,2) £ win(l), k = 2 
and K = (3, 3) € win(2). Then K is compatible with I, for any m2i, 17122 since 
mi3 = —00. Then the interval di5l ) equals [m2i, +oo],/or i = 1 and k = 3. 

Definition 10. An interval relation is an expression Xi G [a, b] + Xk, where a < b 
in T and i,k G [n]. Equivalently, an interval relation is a pair of concatenated 
bivariate linear inequalities Xk + a < Xi < xj. + 6- 

Expressions (ISt- dTTI) imply that the following, at most, four interval relations 
must be true: 



Xk G int(M;i, fc; i, fc) + Xj, i G |/|, A; G |K|. 
Notice that (fT6l ) is tautological when i = k. 



(16) 



Example 3. (Continued) In our running example, take i = 1, k = 2 and I = (1,4). 
The winning pair K = (1,3) is compatible with I and this gives rise to the non- 
empty closed intervals 

int(M;l,2;l,l) = [0,0], int(M; 1, 2; 1,3) = [1,4], 

int(M;l,2;4, 1) = [-5,5], int(M; 1, 2;4,3) = [-4,9] 

and the interval relations 

X3G[1,4]+xi, xi G [-5,5] +X4, X3 G [-4,9] + X4. 

Equivalently, we can write 

xi + 1 < X3 <xi+A, 
X4 — 5 < xi < X4 + 5, 
X4 — 4 < rE3 < X4 + 9. 

Where does, say, the interval relation x^ € [—4, 9] + X4 come from? We know 
that I = (1,4) E win(l) translate into ©-(IZl) and, similarly, K = (1, 3) S win(2) 
translate into 

XI - X3 + 1 = 0, (17) 

-xi + X2 + 1 < 0, (18) 

-xi + X4 - 5 < 0. (19) 

These six expressions imply the concatenated inequalities —4 + X4 < X3 < 9 + X4. 

C 1 

This is possible by compatibility ofK with I, since the tropical minor 



1 



13 attains its value at the main diagonal. 



trop 



Four tropical minors of the maximum matrix M must be checked out, in order 
to decide compatibility of K with /. We can forget repeated minors, keeping just 
one of them. Any minor with repeated columns will be called trivial. Trivial minors 
will be disregarded; they will play no role. A minor is tropically singular if it 
attains its value at both diagonals; otherwise the minor is tropically regular. Of 
course, M{i,k;i,i) is tropically singular. Now, if i ^ k, the minor M{i,k;i,k) 
is tropically singular if and only if the interval relation (fT6l ) reduces to the non- 
tautological bivariate equation 

Xfc = dif (M; i, k)k + Xi_ (= dif (M; i, k)i + Xi). (20) 

Summing up, ii 1 < i < k < m, then the conditions |/| = {ii,i2}' ^ ^ win(i), 
1^1 = {^1)^2}' ^ ^ win(A;) and K compatible with / provide, at most, two 

10 



bivariate linear equations and, at most, four interval relations, namely 

Xi^ = dif(M;ii,i2)i + Xi^, (21) 

Xk^ = dif{M;ki,k2)k + Xk^, (22) 

Xk^ e'mt{M;i,k;ii,ki) + Xi^, (23) 

Xk^ emt{M;i,k-ii,k2) +Xi^, (24) 

Xk^ e int(M; i, k;i2,ki) + Xi^, (25) 

2:^2 G int(M; i, A:; i2> ^2) + a^ia ■ (26) 

The expressions (l23l) and (|25] ) can be combined into one such, using (|2TI ). Sim- 
ilarly, (l24l) and (l26l) can be combined, using (l22l) . In conclusion, expressions (|2T]) - 
(l26ll are reduced to (l2T|>, (EZl) and 

a^fcj G [aM +Xi^, (27) 

a^fcj G [c, d] + Xi^ , (28) 

where a,b,c,d G T, a < b, c < d depend on M, i, k, I and K. In addition, a 
solution X to ([T]) arising from / and from K must satisfy the following: 

Xj < iVi - rriij) ®' (Vk - rukj), (29) 

for all j ^ \I\ U \K\ U F{i) U F{k). This follows from ©. 
Recall that m is the number of rows of the matrix M. 

Definition 11. Let T = (Ji, . . . , /„) ^^ «« m— tuple with I^ G win(/i), /or every 
h G [m]. We i'flj that I is a win sequence /or ([7]) /fJ/i /s compatible with Ii,for all 
1 < i < h < m. 

For a win sequence T = (/i, . . . , /,„), write 

m 

|T|=Ul^^.|- 

Given i, j G |T|, write i ~ j if there exist /c, / G [m] such that i G j/fc|, j G |//| and 
|/fc| n \Ii\ / 0. Closing up under transitivity, we obtain an equivalence relation on 

|T|. 

Definition 12. Let T = {Ii, . . . ,Im) be a win sequence. 
L An index i G [n] is free in Tifi ^ |T|. 
2. An equivalence class for the relation above is called a cycle in T. 



11 



Consider a win sequence T = (Ii, . . . , 1^)- Let c be the number of cycles in 
T. We have 1 < c < card|T| < min{2?7i,n}. After relabeling columns, we can 
suppose that the cycles in T are 

Ci = [hi 

C2 = [k2]\[kil 

Cc = [k,]\%_il 
and that [n] \ [kc] are the free indices, for some 1 < ki < ■ ■ ■ < kc = card |T| < n. 

Theorem 13. Each win sequence T = (Ji, . . . , J^) provides a convex set, solx C 
T", of solutions to the system ([7]). The set solx consists of all the solutions x G T" 
arising from I^^, for all h € [m]. Moreover, 

dim (solx) < n — card |T| + c. (30) 

All solutions to ([7]) are obtained this way. 

Proof. The last statement follows from remark [6l Convexity of solx is trivial, be- 
cause solx is the set of x G T*^ which satisfy a certain system of classical linear 
equations and of inequalities: these come up from ^, for I = Ifi, h £ [m]. Of 
course, the dimension of solx is the dimension of its linear hull. 

For the bound on the dimension, with the notations above, notice that 

• Each cycle C provides a system Tc of bivariate linear inequalities (these are 
related to certain interval relations) only involving Xj, with j £ C (see ex- 
ample 3 in p. O. This can make the dimension decrease (when an interval 
collapses to a point) or stay the same. 

• Each cycle C provides a system Sc of bivariate linear equations only involv- 
ing Xj, with j £ C (see example 3 in p. [141) . There are two possible cases: 
either the only solution to Sc is trivial (i.e., Sc is incompatible over M; see 
example [14] ) or the solution set to Sc is one-dimensional. Notice that, in the 
former case, the condition Xj = — oo, for some j £ C (together with certain 
additional inequalities) might imply x; = — oo for some free index I (this 
explains the expression "at most", in the next item). 

• Each free index / makes the dimension increase by one unit, at most. Notice 
that relations such as ( |29l ) can make the dimension decrease (when yi = —oo 
or yk = — oo) or stay the same. 

There are c cycles and n — card |T| free indices, so that formula (l30l) follows. D 
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Example 14. The set solx can be trivial. Indeed, for 



3 7 


-1 


— oo 








— oo 


— oo 


— oo 


8 


6 7 


— oo 


— oo 




B 


= 


— oo 


— oo 


5 


1 


9 





— oo 








-9 





— oo 


-4 








r 


3 


7 


-1 


8 " 










M 


= 




6 
-9 


7 



5 



1 

-4 









we get 



Then T = ((1, 4), (1, 3), (3, 4)) is a win sequence which implies 

xi — X4 — 5 = 0, 
a;i - X3 + 1 = 0, 
X3 - X4 + 4 = 0, 

but this linear system is incompatible over M. Over T, the only solution is xi = 
x^ = X4 = —oo. In addition, X2 = — 00 because, from the first row of A we know 
that a;2 + 7 < xi + 3, (i.e., X2 — xi + 4 < 0, in normal form); see remove and 
enlarge in p. |6] Here, VL = [4]. 

If no win sequences exist, then the only solution to the system ([T) is trivial. 
The number p of win sequences is no bigger that r*", where r = max{ [^] [^J , n}. 
Even if some win sequence does exist, it may happen that the only solution to the 
system ^ is trivial. 

Example 3. (Continued) Let us finish our current example. The win sequences are 

Ti = ((l,4),(l,3),(3,3))a«JT2 = ((2,4),(l,3),(3,3)). 
For Ti we must solve the systems 



S = Sr,:{ 



Xl 


+ 3 = 


X4 


+ 8, 


Xl 


+ 6 = 


a^s 


+ 5, 


.X3 


+ 1 = 


xz 


+ 1, 



X2 + 7 < Xl + 3, 
X3 - 1 < Xl + 3, 
X2 + 7 < Xl + 6, 
Txi : { X4 + 1 < Xl + 6, 
Xl + 1 < X3 + 1, 
X2 + < X3 + 1, 
^ X4 + 2 < X3 + 1. 
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Writing in normal form, the system S leads to the coefficients matrix 



C 




-1 



(31) 



By Gaussian elimination (which dates back, at least, to the mathematicians of an- 
cient China, in the 2nd or 3rd century B.C., and should better be called the Chinese 
elimination; see /SJ/P- 219), we transform the former matrix into the following up- 
per triangular matrix (also denoted by C) 



C 



1 




(32) 



This gives the following partial result: x^ = X4 + 6 and xi = x^ — 1. Substituting 
X3 by its value (this is usually called backward substitution), we get 



X2, X4 G T, X3 = X4 + 6, xi = X4 + 5. 



(33) 



In particular, the indices 1, 3, 4 belong to the same cycle. 

Writing in normal form and removing tautologies, the system T yields the fol- 
lowing coefficient matrix 



D 



1 


1 








4 


1 





1 





-2 


1 


1 








1 


1 








1 


-5 





1 


-1 





-1 








-1 


1 


1 



(34) 



We want to simplify the information contained in the matrix D as much as pos- 
sible. In order to do so we will produce two matrices E, N such that 

D[x, 1]^ < ^ E[x, if = and N[x, if < 0, 

2 card rows (£^) + cardrows(A^) < cardrows(D). 

Either matrix E or N could be empty. We will denote N by D. The desired simpli- 
fication occurs whenever cardrows(A^) < cardrows(D). 



In order to simplify, we allow the following row operations on the real matrix 



D: 



1. permute two rows, 

2. replace a row by the sum of that row plus another row, so that the resulting 
row has, at most, three non-zero entries. 
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3. remove a row r, when it is superfluous, in the following sense: ifr = (ri , . . . , r„ , a) 
and s = (ri, . . . , r„,, 6) are rows in D, with Vj G M and a < 6 G M, then r is 
superfluous, 

4. remove two rows in D, when they are opposite, and write any one of them into 
the matrix E. 

Warning: multiplication of a row by a negative number is not allowed here. 

By means of such row operations, we can achieve an "almost" upper triangular 
matrix D (the precise term is sub-special; see definition [75l below ) and a matrix E 
(but E turns out to be empty, in this particular case): 



D 



1 


1 








4 


1 





1 





-2 


1 








1 


-5 





1 


-1 





-1 








-1 


1 


1 



(35) 



Notice that row ( — 1, 1,0, 0, 1) in f l34l ) has been eliminated because row (—1, 1,0, 0, 4) 
makes it superfluous (the meaning is that —xi + X2 + 4 < implies —xi + X2 + 1 < 
Oj. Now, the matrix f |J5t stands for the system 

-xi + X2 + 4 <0, 
—xi + X3 — 2 <0, 
—xi + X4 — 5 <0, 
2:2 - 2:3 - 1 <0, 
—x^ + X4 + 1 <0. 

We substitute X3 and xi by their values, as shown in f lJJD . obtaining a new system 
of inequalities, which we must write in normal form and apply row operations to it 
again. Finally, we obtain: 

X2 — X4 — 1 < 0, 



corresponding to the matrix 

L»= [ 1 -1 -1 ] . 
Summing up, the solutions arising from Ti are given by dJJD and 



(36) 



(37) 



X 



X4 + 5 

X4 + 6 
X4 



S.f. X2 — X4 — 1 < 0. 
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Similarly, we get that the solutions arising from T2 are 



, s.t. 2:3 — 6 < X4 < X3 — 3. 



The former example is rather general: we can apply the same procedure to find 
the solutions solx, for any win sequence T. 

It is now convenient to give names to certain types of real matrices that we will 
deal with. 



' X3-1 ' 


X4 + 1 


X3 


3:4 



Definition 15. Let G = {gij) € Mmx{n+i){ 

1. Let G' G Mmxni^) be obtained from G by deleting the last column. 

2. The matrix G is special if each row of G' is a permutation of the n— vector 

(1,-1,0,... ,0). 

3. The special matrix G is super-special if the first non-zero entry of each row 
is L 

4. The special matrix G is sub-special if 

(a) all rows in G are different and different from rows in —G, 

(b) all rows in G' are different, 

(c) ifrow^G', i) = — row(G', k), for some i < k, then k = i + 1, 

9i,n+l < —gi+l,n+l 

i-1 i-i 

and row(G', i) = (0, . . . , 0, 1, 0, . . . , 0, —1, 0, . . . , 0), for some j, I G 

(d) ifrow{G',i) 7^ — row(G',i + 1), for some i, then min{j : gij 7^ 0} < 
min{i : gi+ij ^ 0}. 

Example 16. Matrices ([32\)—/[35]) are special. The following matrix G is sub— special 
but not super-special 



10-103 

-1 1 0-8 

0-1 1 0-4 

1-10 
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Remark 17. Notice that, if G is a non-empty sub-special matrix, then the set of 
X S T" such that G[x, 1]-^ < is neither trivial nor T". The set of all such x will 
be denoted sole. IfG is empty, we define sole as T". 

Remark 18. • By Gaussian elimination, any special matrix C corresponding 
to the system S of equalities C[x, 1] = 0, can be transformed into a super- 
special upper triangular matrix E, such that 

C[x,lf = Q ^ E[x,lf = {), 

cardrows(ii^) < cardrows(C). 
E might be empty. The matrix E is relabeled as C. 

• By the row operations on p. \T4\ any special matrix D corresponding to a 
system T of inequalities D[x,l]'^ < 0, can be transformed into a sub-special 
upper triangular one N and an additional special matrix E, with 

2 card rows (ii^) + cardrows(iV) < cardrows(D). 

Either matrix N or E might be empty. The matrix E corresponds to a system 
of equalities derived from T. The matrix E can be transformed into a super- 
special one. The matrix N is relabeled as D. This procedure is called sub- 
specialization of D. 



3 The algorithm 

An algorithm to solve ([T|) must find first all win sequences. Then, for each win 
sequence T, the algorithm must write two bivariate systems: one system Sr of 
equations and one system Ty of inequalities. Writing these systems into normal 
form and performing remove and enlarge as in p. |6l the algorithm must compute 
a subset f^x ^ [n] and special matrices Cr,Dr G -^mx(n+i)(^^) such that the 
former systems are equivalent to 

Cr[x, if = 0, Dr[x, if < 0, Xj = -oo, j € ftr- (38) 



The set of solutions to (I38I ). denoted soIt = solc^ n soIdt' must now be computed. 
It might be trivial. The non-trivial solutions to ([T|l is the union of solx, as T runs 
over all win sequences. 

ALGORITHM 

• STEP 1 : compute the matrices A, B and M. Replace A and Bhy A and B. 
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STEP 2: Compute all winning pairs, for all i € [m]. Store them in a tridimen- 
sional array W (r rows, 2 columns, m pages). In page i we store all members 
of win(i). Blanks are padded with zeros. 

STEP 3: Compute all win sequences. Store them in a tridimensional array 
WS {m rows, 2 columns, p pages), with < p < r™. No entry of WS is 
zero. If WS is empty, then the only solution to system ([T) is trivial, RETURN. 

FOR each win sequence T 

- STEP 4: Compute the set fix Q [n] and the special matrices Ct and 
Dr. 

- STEP 5: Sub-specialize matrix Dx- In order to do so, first, obtain 
special matrices E^, Ny such that 

Dr[x, if <0 ^ Er[x, if = and Nr[x, if < 0, 

2cardrows(ii^x) + card rows (A^x) ^ cardrows(L'x)- 

Either matrix Ey or A'^x could be empty. Denote A^x by Dy . By row 
operations, work on D~f to make it sub-special. 

- STEP 6: Concatenate the matrices Cx and Ey into a matrix, which we 
can denote again by Cx- By Gaussian elimination, work on Cx to make 
it super-special and upper triangular. The number of rows of Cx may 
decrease. Solve the classical linear system Cr[x, if = 0. The set J7x 
may enlarge. The solution set in T" is denoted solc^ ; it may depend on 
a number of parameters. If solc-f is trivial, GO TO WORK WITH THE 
NEXT WIN SEQUENCE. 

- STEP 7: If Dy is non-empty, substitute x G solc^ into -Dx [x,lf ^ 
to obtain a new system of linear inequalities. Write this system in nor- 
mal form and denote it Dt[x, if < again. The set ilx may enlarge. 

- STEP 8: If Dr is non-empty, sub-specialize matrix Dy. Obtain special 
matrices E~f, Nj- such that 

Dr[x, if <0 ^ Er[x, if = and Nr[x, if < 0, 

2cardrows(ii^x) + card rows (iVx) < cardrows(L'x)- 

Either matrix Ey or A'^x could be empty. Denote A^x by Dy. By row 
operations, work on D~f to make it sub-special. If Ej- is empty, then 

solx = solcr n solDy ■ 
Otherwise, GO TO STEP 6. 
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• END FOR 

All the solutions to ([B are UxeH'S solx- 

We have programmed the former algorithm to solve system ^. Working over 
Q, let us compute the complexity of it. The arithmetic complexity counts the num- 
ber of arithmetic operations (+, — , max, min, <, = and >, in our situation) in the 
worst possible case. 

Our programme is divided into two parts. In the first part, we determine all the 
win sequences. Say we get p win sequences. The arithmetical complexity of this 
partis 

0{m n p). 

In the second part, we compute the matrices Cr , Dy and all the solutions (if any), 
for each win sequence T. The arithmetic complexity of the second part is 

0{m{m +n)p). 

Since the maximum number of winning pairs is r = max{[^] [§J,n}, then p < 
r™, where r is O(n^). This gives an exponential arithmetical complexity! But, let 
us take a closer look. Clearly, the bigger n, the more winning pairs we may have, for 
each i G [in]. On the other hand, the bigger m, the fewer win sequences we have, 
in probability, due to the compatibility requirement. Indeed, given winning pairs 
/ G win(i), K S win(A;) with 1 < i < A; < rri, let us assume that the probability 
of K being compatible with I as 1/2 (this assumption is rather reasonable, since 
"K is compatible with /" is a yes/no event). Thus, given any sequence of pairs 
T = (/i, . . . , Im), the probability of T being a win sequence is, roughly, 

1 1 



2(2) 2'-^ ■ 

This proves that if m is big, then we expect p rather small. In particular, the worst 
case, p = r"^,is unlikely to happen. With this in mind, an average complexity (see 
II2I ) for the first part is 

0(m2n3+2'»/2'"') = o(m22(3+2™)log2'^-"*2) 

and it will be, at most polynomial 0{m'^), if log2 n < 3^^ ~ y- ^'^^ ^^ second 
part we get two terms: 

0(m3n2"^/2""') = 0(m322'"l°g2"-"^') 

and it will be, at most polynomial 0{m^), if log2 n < y , and 

0(mni+2™/2"^') = 0(m2(^+2'")l°g2'^-™') 

2 

and it will be, at most polynomial 0{m), if log2 n < ^^m ~ T' 
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4 Some examples 

Example 19. Given 

A = 
we get 



The only win sequence is T = ((2, 3), (1, 1), (2, 1)). The solutions arising from T 
are 

X = X3 
X3 

Here we have m = n = r = 3, r"^ = 27, r'"/2"* ~ 0.0527 and p = 1. The 
programme execution lasted 0.038883 seconds with Matlab R2007b in a portable 
PC working with Windows Vista. 

Example 20. (From EJTi) Given 

^ _ _ _ . , 2 6 

^" - " -oo 2 

-oo 
-oo 



B 



we get 



— oo 


- 


— oo 
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oo 
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oo 
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oo 4 


5 




6 — oo — oo 
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— oo — oo — 
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oo 
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" 
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5 4 2 6 


3 


5 


6 2 4 2 



oo 
4 



The win sequences are T I = ((4, 1), (2, 1)), T2 = ((4,3), (2, 1)), T3 = ((5, 1), (2, 1)), 
T4 = ((5, 3), (2, 1)), T5 = ((6, 1), (2, 1)), Te = ((6, 3), (2, 1)), T7 = ((7, 1), (2, 1)) 
an^T8 = ((7,3),(2,l)). 

Here we have m = 2, n = 7, r = 12, r™ = 144, r"'/2'^' = 9 and p = 8. The 
programme execution lasted 1.210653 seconds with Matlab R2007b in a portable 
PC working with Windows Vista. 
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The solutions arising from Ti are 



X4 
X4-2 

3:4 

Xq 
X-j 



s.t. 



X3 — X4 + 5 < 0, 
-X4 + X5 + 4 < 0, 
-X4 + xg + 2 < 0, 
-X4 + rc7 + 6 < 0. 



The solutions arising from T2 are 



X2 + 2 

X4 — 5 
a;4 

xg 



2:2 + 2 < X4 < X2 + 4, 

-X2 + 2:5 - 3 < 0, 

-X2 + 2;6 - 1 < 0, 
5.f. — X2 + Xy — 3 < 0, 

-X4 + X5 + 4 < 0, 
— X4 + xg + 2 < 0, 
— X4 + X7 + 6 < 0. 



The solutions arising from T3 are 

X5 + 4 
X5 + 2 

2:3 

X4 

2:5 

X6 
X7 

r/ie solutions arising from T4 are 

X2 + 2 

X5 - 1 

X6 
X7 



S.t 



2^3 - 2:5 + 1 < 0, 

X4 — X5 — 4 < 0, 
-X5 + xg - 2 < 0, 

-X5 + X7 + 2 < 0. 



X2 - 2 < X5 < X2, 

X4 — X5 — 4 < 0, 
-X2 + X4 — 5 < 0, 

-X2 + X6 - 1 < 0, 

-X2 + X7 — 3 < 0, 
-X5 + X6 - 2 < 0, 

-X5 + X7 + 2 < 0. 
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The solutions arising from T5 are 

xq + 2 

X6 

X3 

X4 , s. t. 

X5 

X6 
Xj 

The solutions arising from Tg are 



a^a - 2;6 + 3 < 0, 

X4 — xq — 2 < 0, 

X5- xq + 2 <0, 

-xq + X7 + 4 < 0. 



X2 + 2 

X2 

xq-3 

X4 , s.t. 

X5 

Xq 
Xj 



X2 < Xq < X2 + 1, 
X4 — Xg — 2 < 0, 

X5- Xq + 2 <0, 
-X2 + 3:4 — 5 < 0, 
-X2 + 2:5 - 3 < 0, 
-X2 + 2:7 — 3 < 0, 

-XQ + Xj + 4:< 0. 



The solutions arising from T7 are 

X7 + 6 

X7 + 4 

X3 

X4 , s.t. 

X5 

X6 
Xj 

The solutions arising from Tg are 



X3 — x-j — 1 < 0, 
X4 — 3:7 — 6 < 0, 
X5 — XT — 2 < 0, 
xg — X7 — 4 < 0. 



X2 + 2 

X2 
X7 + I 

X4 , S.t. 

X5 

Xq 
Xj 



X2- 4:< Xj < X2-2, 
X4 — X7 — 6 < 0, 

X5 — X7 — 2 < 0, 

xg — X7 — 4 < 0, 
-X2 + X4 — 5 < 0, 
-X2 + X5 - 3 < 0, 

-X2 + Xg - 1 < 0. 



The set of solutions has pure dimension 5. 
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